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vjded by cosmic rays and radioactive nuclei. For the purpose of our 
calculations we assume a solar protoplanetary disk to be an a disk 
with the opacity law taken from Ruden and Pollack [3] and the 
ion\zation state taken from Stepinski [4]. We choose a - 0.08, and 
M- 1(H M© per year. We assume a disk surrounding a 1-M© star 
and extending from 0.2 AU up to 40 AU. Figure 2 shows the time 
evolution of the magnetic Field from the initial Field Q(r) -0.1 in 
units of the equipartition value at r - 40 AU. At First the Field 
increases. sharply at the inner radii, decays at the middle radii, and 
remains unchanged at the outer radii. By the time t « 10 yr, the 
magnetic field in the innermost portion of the disk achieves equilib- 
rium. As time progresses the magnetic Field achieves equilibrium at 
larger and larger portions of the inner disk. At the same time, the 
field continues to decay at the middle radii, but the decaying region 
shifts outward as a result of radial diffusion, and the magnetic Field 
in the outer parts starts to show some growth. By the time t - 1 00 y r 
the whole region within 3 AU has reached equilibrium. Radial 
diffusion from the regions of strong magnetic field stops the further 
decay of the field within the region where the local growth rate is 
negative, and the Field is now actually growing there. The magnetic 
field in the outer parts of the disk continues to grow. By the time t - 
2000 yr, the magnetic field in almost the entire disk has reached 
equilibrium. Total equilibrium is achieved at roughly t - 4400 yr. 
The final configuration of the magnetic Field follows closely the 
distribution of the equipartition value magnetic field, except at the 
middle radii. 

Conclusions: The fipal configuration of a dynamo-generated 

magnetic field is independent of unknown initial conditions. How- 
ever, initial conditions influence the way the magnetic field evolves 
toward its equilibrium, as well as the time needed to achieve such 
equilibrium. Evolution from initial conditions without field rever- 
sals (presented here) leads to an equilibrium field in a time that is 
very short in comparison with^disk viscous time. Evolution from 
initial conditions with field reversals (not shown here) leads to an 
equilibrium in a time 10— 10 2 times longer, as radial diffusion 
destroys field reversals. In equilibrium, the field has a magnitude of 
the order of the equipartition with the kinetic energy of turbulence. 




Fig. 2. Time evolution of the magnetic field in a proto planetary disk from the 
initial condition Q - 0.1 at t - 0 represented by the horizontal solid line. 
Magnetic field Q is measured in un i ts of B 0 - ( 40 AU ) . Radial di stance from 

the central star is measured in AU. 


Such a Field could have a substantial effect on the structure and 
dynamical evolution of thin disks. From an observational point of 
view, the magnetic Field is concentrated close to the inner disk’s 
radius, so it could be difficult to distinguish it from a stellar mag- 
netic Field, provided that a central star has a strong field. 
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Fig. 1. Radial distribution of magnetic field Q is plotted against dimensionless 
radius r/r^ at various times for the case of an accretion disk around a compact 
star. The plots (i-f), in order of increasing time, correspond to t - 1 0, 1 0 2 , 1 0 3 , 
10* 10 5 , lO^and l0 7 s respectively. The dottedline shows the radial distribution 
ofB^. After about I - 1 0 7 s the magnetic Field equilibrates everywhere in a disk 
at about half the equipartition value. \ 
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We consider the structure and emission spectra of nonthermal 
accretion disks around both strongly and weakly magnetized neu- 
tron stars. Such disks may be dissipating their gravitational binding 
energy and transferring their angular momentum via semicontinuous 
| magnetic reconnections. We consider specifically the structure of 
l the disk-stellar magnetospheric boundary where magnetic pressure 
; | balances the disk pressure. We consider energy dissipation via 
* 1 reconnection of the stellar field and small-scale disk turbulent Fields 
e of opposite polarity. Constraints on the disk emission spectrum are 
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c i The nonaxisymmetric stability of self-gravitating, geometrically 

/thick accretion disks has been studied for protostellar systems 
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having a wide range of disk-to-central object mass ratios. Global 
eigenmodes with four distinctly different characters have been 
identified using numerical, nonlinear hydrodynamic techniques. 
The mode that appears most likely to arise in normal star formation 
settings, however, resembles the “eccentric instability” that has 
been identified earlier in thin, nearly Keplerian disks: It presents an 
open, one-armed spiral pattern that sweeps continuously in a trail- 
■ ing direction through more than 2 -n radians, smoothly connecting 
the inner and outer edges of the disk, and requires cooperative 
motion of the point mass for effective amplification. This particular 
; instability promotes the development of a single, self- gravitating 
* clump of material in orbit about the point mass, so its routine 
appearance in our simulations supports the conjecture that the 
1 eccentric instability provides a primary route to the formation of 
"short-period binaries in protostellar systems. 
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\ We have developed a tool to solve the radiative transfer equation 

‘for a three-dimensional astrophysical object on the SIMD computer 
MasPar MP- 1 . With this tool we can rapidly calculate the image of 
;such an object as seen from an arbitrary direction and at an arbitrary 
: wavelength. Such images and spectra can then be used to directly 
'compare observations with the model. This tool can be applied to 
" many different areas in astrophysics, e.g., HI disks of galaxies and 
" polarizedradiative transfer of accretion columns onto white dwarfs. 

Here we use this tool to calculate the image and spectrum of a simple 
~ model of an accretion disk. 
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j The study of magnetized plasmas in astrophysics is complicated 

l by a number of factors, not the least of which is that in considering 
? magnetic fields in stars or accretion disks, we are considering 
' plasmas with densities well above those we can study in the labo- 
* ratory. In particular, whereas laboratory plasmas are dominated by 
"the confining magnetic field pressure, stars, and probably accretion 
disks, have magnetic fields whose P (ratio of gas pressure to mag- 
netic field pressure) is much greater than 1 . Observations of the Sun 
Suggest that under such circumstances the magnetic field breaks 
apart into discrete flux tubes with a small filling factor. On the other 
hand, theoretical treatments of MHD turbulence in high-p plasmas 
lend to assume that the field is more or less homogeneously distrib- 
uted throughout the plasma [ij. 

' Here we consider a simple model for the distribution of magnetic 
flux tubes in a turbulent medium. We discuss the mechanism by 


which small inhomogeneities evolve into discrete flux tubes and the 
7 size and distribution of suchHux tubes. We then apply the model to 
; accretion disks. We find that the fibrilation of the magnetic field 
' does not enhance magnetic buoyancy. We also note that the evolu- 
~ tion of an initially diffuse field in a turbulent medium, e.g., any 
* uniform field in a shearing flow, will initially show exponential 
growth as the flux tubes form. This growth saturates when the flux 
tube formation is complete and cannot be used as the basis for a self- 
sustaining dynamo effect. Since the typical state of the magnetic 
field is a collection of intense flux tubes, this effect is of limited 
interest. However, it may be important early in the evolution of the 
galactic magnetic field, and it will play a large role in numerical 
simulations. Finally, we note that the formation of flux tubes is an 
~ essential ingredient in any successful dynamo model for stars or 
accretion disks. 

We will consider an idealized situation in which there exists a 
turbulent cascade with a scale L and a turbulent velocity, on the 
scale of V T . We will assume that the magnetic field has an rms 
Alfvfch speed where ~ V T . We wi II also assume that the typical 
scale of curvature for the field lines is L. These assumptions are less 
restrictive than they may appear. If the turbulent cascade actually 
extends to larger length scales and higher velocities, then the mag- 
netic field is dynamically insignificant on these larger scales and we 
can still confine, our attention to scales of size L or smaller. If the 
magnetic field is in a shearing flow, surrounded by turbulence of its 
own creation, thenthe near equality of V T and V A is guaranteed, as 
well as the curvature of the magnetic field lines on the scale L. 

The field lines wilLtend to stretch at a rate -Vj/L. If the plasma 
is highly conducting then the same amount of matter will be en- 
trained on a progressively iQngerand longer flux tube. In a stationary 
state this stretching will be balanced by the pinching off of closed 
loops. These loops will have a radius ~L and a longitudinal compres- 
sive force - pV^/L. This tension will be opposed, usually by turbu- 
lent stretching with a force of ~V£/L, Some large fraction of the time 
the loops will collapse. Regardless whether the internal pressure of 
the loop is dominated by the magnetic field or gas pressure the 
magnetic tension will decrease more slowly than the turbulent 
stretching force and the loop will collapse to a plasmoid ball, whose 
energy is slowly lost to microscopic dissipation. This process will 
tend to remove matter from the flux tubes at a rate of -VyL, which 
is rapid and will produce largely evacuated flux tubes under almost 
any circumstances. If we start from a uniform or nearly uniform 
field, this process will end when the same amount of flux is divided 
into some number of intense flux tubes with a magnetic pressure 
equal to the ambient pressure and a local P of order unity or less. The 
final mis Alfvdn velocity will be the geometric mean between its 
initial value and the local sound speed. This increase will occur at 
a rate comparable to V-f/L, in agreement with the results of numeri- 
cal experiments [2,3]. 

What will be the typical radius of the individual flux tubes? A 
single flux tube with an internal AIfv£n speed of V A , ~ c s . and 
exposed to an ambient turbulent velocity of V T , will remain coupled 
to the fluid provided that r, < (V 1 7V At ) 2 L. On the other hand, these 
tubes will impede the flow, and thereby reduce the ambient fluid 
velocity below V T , if the total number N is large enough that Nr t /L 
Is greater than 1. The requirement that the magnetic energy be 
divided into N flux tubes is just the requirement that Nr j V V £L 2 , 
which implies that the flux tubes will not impede the flow if r t is 
comparable to, or greater than, L(V A /V At ) 2 . We conclude that the 


